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Abstract. We consider the equilibrium equations for a conducting elastic rod placed 
in a uniform magnetic field, motivated by the problem of electrodynamic space tethers. 
When expressed in body coordinates the equations are found to sit in a hierarchy of non- 
canonical Hamiltonian systems involving an increasing number of vector fields. These 
systems, which include the classical Euler and Kirchhoff rods, are shown to be completely 
integrable in the case of a transversely isotropic rod; they are in fact generated by a Lax 
pair. For the magnetic rod this gives a physical interpretation to a previously proposed 
abstract nine-dimensional integrable system. We use the conserved quantities to reduce 
the equations to a four-dimensional canonical Hamiltonian system, allowing the geometry 
of the phase space to be investigated through Poincare sections. In the special case 
where the force in the rod is aligned with the magnetic field the system turns out to 
be superintegrable, meaning that the phase space breaks down completely into periodic 
orbits, corresponding to straight twisted rods. 

AMS classification scheme numbers: 74K10, 78A30, 70H06 
1. Introduction 

The problem of an elastic conducting wire in a magnetic field is a classical one in 
magnetoelasticity. For instance, it is well-known that a straight current-carrying wire 
held between pole faces of a magnet buckles into a coiled configuration at a critical 
current [ I , Section 10.4.3]. A rigorous bifurcation analysis of this buckling problem (for 
a uniform magnetic field directed parallel to the undeformed wire) was developed in a 
series of papers by Wolfe. He first considered a nonlinearly-elastic string model for the 
wire, i.e., a perfectly flexible elastic line, and found that an infinite number of solution 
branches bifurcate from the trivial straight solution [2, 3], much like in the Euler elastica 
under compressive load. The equations can in fact be solved exactly and the non-trivial 
solutions are exact helices. In subsequent work Wolfe modelled the wire as a rod [4, 5, 6]. 
In addition to extension a rod can undergo flexure, torsion and shear, and for the case of 
welded boundary conditions it was found that in certain cases bifurcation occurs, again 
with an infinity of non-trivial equilibrium states. In all these studies the wire is assumed 
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to be non-magnetic. Magnetoelastic instability of a current-carrying rod due to its own 
magnetic field has also been considered [7, 8, 9]. 

The problem of a conducting rod in a magnetic field has recently attracted interest as 
a model for electrodynamic space tethers, i.e., conducting cables that exploit the earth's 
magnetic field to generate thrust and drag (Lorentz) forces for manoeuvring [10, 11]. 

Here we consider the hitherto unexplored Hamiltonian structure and integrability of 
the equilibrium equations for a rod in an external magnetic field. Like Wolfe we assume 
the magnetic field to be uniform and we shall concentrate on the usual case of a rod 
that is uniform, intrinsically straight, inextensible and unshearable. For this case the 
equilibrium equations formulated in body coordinates are found to sit in a hierarchy of 
rod models, described by an increasing number of three-dimensional field vectors, with an 
interesting non-canonical Hamiltonian structure. The first two members of this hierarchy 
are the classical single-field force-free (Euler) rod and the two-field Kirchhoff rod (the 
two fields being the force and moment in the rod). Both these rods have well-known 
dynamical analogues in, respectively, the swinging planar pendulum and the spinning top 
[12, 13]. The 'magnetic rod', the next level in the hierarchy, does not seem to have a 
dynamical analogue. However, it does give a physical realisation to (a special case of) the 
'twisted top' introduced in [11]. The equations of this top were introduced as an abstract 
mathematical construction obtained by extending the Lie-Poisson bracket of the heavy 
top by a cocycle (Leibniz extension). The resulting three-field bracket is the simplest non- 
semidirect extension of the heavy-top bracket [15] and the authors in [14] were interested 
in how this extension affects the dynamics of the system. We show that the 'twisted top' 
is in fact a magnetic rod. 

By using the Casimir invariants for the magnetic rod we perform an explicit reduction 
to a six-dimensional canonical Hamiltonian system with a two first integrals in addition 
to the Hamiltonian. This reduction is valid for all but an invariant subspace of the phase 
space where the magnetic field is aligned with the force in the rod. 

The magnetic rod can be regarded as a perturbation of the Kirchhoff rod and is found 
to retain the latter's property that the equations are completely integrable if the rod is 
transversely isotropic, i.e., if the principal bending stiffnesses are equal (the dynamical 
analogue is the symmetric Lagrange top). In fact, the subhierarchy of isotropic rods is 
generated by a Lax pair. On the alignment subspace the isotropic rod is found to be 
superintegrable and the equations can be solved directly by algebraic means. All solutions 
correspond to twisted straight rods. 

The next (four-field) level of the hierarchy also has physical meaning: it describes a 
rod in a special type of non-uniform magnetic field generated by what might be called a 
uniform hypermagnetic field. Alignment of these two fields again defines a superintegrable 
case with a breakdown of the phase space into periodic orbits corresponding to twisted 
straight rods. The mathematical structure of the hierarchy extends beyond this four-field 
model but the physical meaning becomes less clear. 

The outline of the paper is as follows. In Section 2 we set up the geometrically exact 



Integrability of a conducting elastic rod in a magnetic field 



3 



Cosserat theory used for modelling a rod. The hierarchy, with increasing number of fields, 
is presented in Section 3. We briefly review the classical cases in the present context and 
introduce the magnetic rod in Subsection 3.3. We show that the isotropic case is integrable 
in the sense of Arnold-Liouville. In Section 4 we perform the reduction of the magnetic 
rod equations to a six- dimensional canonical Hamiltonian system. We also identify the 
superintegrable case. Section 5 presents the hypermagnetic rod, while Section 6 gives the 
Lax pair for the subhierarchy of isotropic rods. 



2. Cosserat rods 



Cosserat theory for slender rods describes a rod as an arbitrary curve in space endowed 
with mechanical properties such as bending and torsional stiffnesses [ ] ]. Forces and 
moments, whose balance equations will be introduced in the following sections, are 
assumed to be suitable averages over the rod's cross-section and to act at the centreline. 
The configuration of a rod is defined by a vector function r (s), describing the centreline of 
the rod, and an attached right-handed orthonormal triad of directors {di(s), d 2 (s), d 3 (s)} 
describing the orientation of the cross-section. We assume the rod to be inextensible and 
unshearable. In mathematical terms this is expressed by the conditions \r'\ = 1 and 
r' • d\ = r' • di = 0, respectively, where ()' = d/ds. Thus we can write r' = d% and the 
parameter s measures arclength along the rod in any configuration. It will be convenient 
to express components of vectors with respect to the director (or body) frame; for any 
vector v the triple of components Vi — v ■ di (i — 1, 2, 3) will be denoted by the sans-serif 
symbol v. 

On introducing a right-handed orthonormal frame {ei,e2,e^} fixed in space we can 
write 

di = Rei (i= 1,2,3), (2.1) 

where R is a rotation matrix. Differentiation gives 

d\ = R' ei = R'R~ l di =:ux d t , (2.2) 

where u is the vector of generalised strains. The body components «j = u ■ di are the 
curvatures (i = 1, 2) and the twist {i = 3) of the rod, for which, from (2.2), we can write 

= \ e ijk d 'j ■ (2.3) 

where e^k is the standard Levi-Civita permutation symbol. 

Constitutive relations specify how the strains relate to the stresses on the rod. Here we 
assume the rod to be uniformly hyperelastic, i.e., we assume that there is a strain energy 
function W(u) and that the stresses (body moments) mj can be obtained by partial 
differentiation: 

dW 

rrii = 

OUi 
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We shall only consider the quadratic form 

W (u) = \k x u\ + \k 2 u\ + \k,u\ (2.4) 

for a uniform, intrinsically straight and prismatic rod, corresponding to the linear 
constitutive relations (Hooke's law) 

mi = Kiin, m 2 = K 2 u 2 , m 3 = K 3 Us, (2.5) 

where K\ and K 2 are the bending stiffnesses, about d\ and d 2 respectively, and K% is 
the torsional stiffness, about the tangent d 3 . The rod is called transversely isotropic (or 
isotropic for short) if K\ = K 2 =: K. For rods of homogeneous cross-section we have the 
relations: K\ = Eli, K 2 = EI 2 , K 3 = GJ, where E is Young's modulus, G the shear 
modulus, I\ and I 2 the second moments of area of the cross-section about d 2 and d±, 
respectively, and J the polar second moment of area. 

3. A hierarchy of rod models 

3.1. The force-free rod - Single-field model 

The equilibrium equation for a force-free rod is 

m! = 0. (3.1) 
In the director frame the equation can be written as a non-canonical Hamiltonian system 

m' = J(m)VH, (3.2) 

where the skew-symmetric structure matrix J = J (m) is given by 

' — m3 m 2 ' 
J = —J T = m 3 —mi ='■ rh (3.3) 
^ — m 2 mi j 

and the Hamiltonian is 

ft = ^m-u, (3.4) 

with the Ui as given in (2.5). For any two functions of m we can also introduce the 
Lie-Poisson bracket 

{/, 0}(m) = - m ■ (V m / x V m g) , (3.5) 
twist 

and write (3.2) as 

m' = {m,ft} (m) = m x u, (3.6) 

which, on using (2.5), is also Euler's equation for the rotational motion of a free rigid 
body about its centre of mass. 

For a general Hamiltonian system with Poisson bracket { • , • } and Hamiltonian TC a first 
integral is a function I satisfying {I, Ti} = 0. Two integrals I\ and I 2 are said to be in 
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involution (or to commute) if 1%} = 0. In particular, if the bracket is non-canonical 
there exist conserved quantities, called Casimirs, that are in involution with any smooth 
function. Equivalently, C is a Casimir if VC is in the null-space of the structure matrix 
J [17]. An m-dimensional non-canonical Hamiltonian system is completely integrable in 
the sense of Arnold-Liouville [18] if it possesses k Casimirs and I first integrals (including 
the Hamiltonian), all in involution and with linearly independent gradients, such that 

m = 21 + k. (3.7) 

In particular, the solution of completely integrable systems can be reduced to quadrature. 
The null-space of the structure matrix (3.3) is spanned by 

1 T 

VCi = - (m 1 ,m 2 ,m 3 ) 

and hence 

Ci = m • m (3.8) 

is a Casimir. It describes the fact that the magnitude of the total moment is constant along 
the rod. Since TC is an integral, it follows from (3.7) that (3.6) is completely integrable, a 
well-known fact from rigid-body dynamics. If G\ = then the system is degenerate. 

The force-free rod is not merely integrable, it is superintegrable, meaning that it has 
more independent integrals, r, than the number of degrees of freedom, d. Although 
independent, the integrals will not be in involution. In general, in a superintegrable 
system 2d — r conserved quantities define r-tori in 2<i-dimensional phase space. A system 
is minimally superintegrable if solutions exist on (d— l)-tori and maximally superintegrable 
if solutions exist on 1-tori [19]. The force-free rod has three degrees of freedom and 
there are four independent integrals: the Hamiltonian 7i and the three components 
of the moment m in the fixed frame (from (3.1)). Hence, phase-space trajectories are 
restricted to invariant 2-tori within the Liouville 3-tori [, )] and the system is minimally 
superintegrable. In the isotropic case k 2 = (m\ + rn%)/K 2 and r = m^/K^ are constant 
and satisfy K 2 k 2 + K 2 t 2 = Ci, and the rod configurations are twisted helices of curvature 
n and twist r [21]. 

In the special case where K\ = K 2 = K3 = K the equation in the body frame 
becomes m' = 0. This introduces a further, fifth, independent integral (we can choose the 
Hamiltonian, two components of the moment m in the fixed frame and two components 
of m in the body frame as independent integrals) and the phase space is foliated by 
1-tori. This case is maximally superintegrable. The corresponding rod configurations are 
resonant helices that have the material frame moving with the Frenet frame, with the 
twistless ring (r = 0) and the twisted straight rod (k = 0) as extreme cases. 

3.2. The Kirchhoff rod - Two-field model 



The equilibrium equations are 

m! + r' x n = 0, n' = 0, 



(3.9) 
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where m and n are the moment and force in the rod. In the director frame the equations 
can be written in Hamiltonian form 




= J (m, n) VH (m, n) , (3.10) 
where the structure matrix J = J (m, n) is given by 

'-^-(::) en) 

and the Hamiltonian is now 

H = ■ m + d 3 • n, where d 3 = (0,0,1). (3.12) 

Canonical and non-canonical Hamiltonian formulations of this problem have also been 
discussed in [13]. Alternatively, we can write 

m' = {m, H} {m n) = m x u + n x d 3 

n ' = { n ^}(m,n) = n x u, 
in terms of the Lie-Poisson bracket on (m, n) given by 

{/, <7}(m,n) = - m • (Vm/ x V m g) - n • (V m / x V n g + VJ x V m g), (3.14) 

s -v ' 

force 

where an extra (semidirect) term has been added compared to (3.5). The equations (3.13) 
are also those for the motion of a heavy top. 

The null-space of the structure matrix J is two-dimensional and spanned by 





VCi = and VC 2 = - (3.15) 



and hence the Casimirs are 

d = n • m, (3.16) 
C 2 = n • n. (3.17) 

The Casimir C\ describes the conservation of the moment about the force vector, while 
C*2 describes the conservation of the magnitude of force in the rod. If C 2 — then the 
rank of the structure matrix changes and we recover the previous case of the force-free 
rod. 

Casimirs are conserved quantities that only depend on the structure of the balance 
equations; first integrals are dependent on the specific constitutive relations chosen. In 
addition to the Hamiltonian and the two Casimirs a first integral is required if the system 
is to be completely integrable. There are two cases, both well-documented: 

The Lagrange case has an integral given by 

h = Km-d 3 (if K 1 = K 2 = K), 

Thus, if the two bending stiffnesses are equal then the twist m 3 is a conserved quantity. 
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The Kovalevskaya case has an integral given by 

h = (K\m\ - Klm\ + n 3 f + {2K 1 K 3 m 1 m 3 - mf (if K x = K 3 = 2K 2 ). (3.18) 

The above condition on the bending stiffnesses renders the Kovalevskaya rod 
somewhat unphysical since it correspomds to a negative Poisson ratio. This ratio, 
defined as v = E/2G — 1, is positive for ordinary materials. However, novel materials 
with negative effective Poisson ratio are now well known. For instance, experimental 
measurements of bending and torsional stiffnesses of DNA molecules have led to the 
generally accepted range 0.7 < K/K 3 < 1.5 [22]. 

There is another case which is not completely integrable on the entire phase space but is 
dependent on the value of Casimir C\. 

The Chaplygin-Goryachev case requires that the initial conditions must satisfy 
m • n = 0; 

then we have the integral 

h = K 2 m 2 {Klm\ + K 2 2 m 2 2 ) - K 3 m 3 n 2 (if K t = AK 2 = K 3 ). 

It is simple to show that the integrals are in involution with respect to the bracket in 
(3.14) and that 7i = h, C\ = a and I\ = (3 define three-tori on which configurations 
generically exist. 



3. 3. The conducting rod in a uniform magnetic field - Three-field model 



Now consider a rod placed in a uniform magnetic field B. The rod carries a uniform 
current J = Ir' of strength / along the centreline (we assume the rod to be sufficiently 
slender for eddy currents within the cross-section to be ignorable). The rod then 
experiences a Lorentz body force Fl = I x B = Ir' x B = Id 3 x B and the force 
balance equation becomes n' + F ^ = 0. We shall write IB = B so that the equilibrium 
equations take the symmetric form 

n' + r' x B = 0, B' 



m 



r x n 



0. 



0. 



(3.19) 



In writing down these equations we have assumed that the current in the rod is moderate 
so that the effect of the magnetic field generated by the current is negligible compared 
to the external magnetic field. Eqs (3.19) have been considered by Wolfe in bifurcation 
studies from the trivial straight solution [4, 5, 6]. 

In the director frame the governing equation is a non-canonical Hamiltonian system of 
the form 

/m V 



n 
B 



J (m, n, B) VH (m, n) 



(3.20) 



V b / 
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where the structure matrix J 



J (m, n, B) 



is given by 





f 


m 


n 


B\ 


J = -J T = 




n 


B 







\ 


B 





0/ 



and the Hamiltonian is 
1 



n 



m ■ u + dc 



(3.21) 



(3.22) 



We note that the Hamiltonian is the same as for the Kirchhoff rod in (3.12): the effect of 
the magnetic field is only present in the structure matrix. The governing equations can 
be written as 



m' = {m, W} (m n B) = m x u + n x d 3 , 
n' = {n,W} (m n B) = n x u + B x d 3 , 
B' = {B,%, n , B) = Bxu ) 
where the Lie-Poisson bracket on (m, n, B), 

{/^}(m,n,B) = ~ m • ( V m/ X V m #) - n ■ (V m / X V n £ + V„/ X V m #) 

- B • (V m / x V B g + V B / x V m g) - B ■ (V n / x V n < 



(3.23) 
(3.24) 
(3.25) 



evolution of field effect of field 

has been extended from the previous bracket (3.14) by the addition of two more terms. 
The first term, due to the evolution of the magnetic field in the director frame, does not 
affect the governing equations since the Hamiltonian is independent of B. The second term 
contains all the magnetic effects. This term makes the bracket extension non-semidirect 
[15]. Note that (3.25) is merely the trivial equation B' = rewritten in body components. 
Eqs (3.23)-(3.25) give a physical realisation of the 'twisted top' introduced as an abstract 
mathematical system in [14]. 

There are three Casimirs, given by 

d = ■ n + m • B, (3.26) 

C 2 = B • n, (3.27) 

C 3 = B • B. (3.28) 

We observe that the magnitude of force is no longer conserved, but as a result of rotational 
symmetry, the force component in the direction of the magnetic field is conserved resulting 
in (3.27). The magnitude of the magnetic interaction is conserved, thus (3.28). Casimir 
(3.26) however does not seem to have a physical interpretation. 
In the isotropic case (K\ = K 2 ) we have the integrals 

ii = firm-d.3 (if Ki = K 2 = K), (3.29) 
J 2 = n ■ m + KB ■ d 3 (if K x = K 2 = K). (3.30) 

The first of these expresses conservation of twist in the rod, as in the Kirchhoff case. The 
second integral, like does not seem to have a physical interpretation when C 3 7^ 0. We 



Integrability of a conducting elastic rod in a magnetic field 



9 



note that the Lagrange integrability condition K\ = K 2 is unaltered by the magnetic 
field. The authors in [14] give numerical evidence (in the form of chaotic orbits) that the 
same is not true for the Kovalevskaya case: the magnetic rod with K\ = K3 = 2K 2 is not 
integrable. Of course, a .B-perturbed condition on the stiffnesses may exist for which the 
system is integrable. 

It is a straightforward task to check that all the integrals (3.22), (3.29) and (3.30) are 
independent and in involution with respect to the Lie-Poisson bracket (3.26). Hence in 
the isotropic case the system is completely integrable in the sense of Arnold-Liouville. In 
the next section we use the conserved quantities to reduce the nine-dimensional system 
(3.23)-(3.25). 



4. Reduction of the magnetic rod 

In this section we use the three Casimirs (3.26)-(3.28) to reduce the nine-dimensional 
non-canonical Hamiltonian system to a six- dimensional canonical Hamiltonian system. 
This is possible (at least locally) provided the structure matrix J is of constant rank 
everywhere [17, Section 6.2]. We perform the reduction by constructing a coordinate 
transformation from the nine coordinates (m, n, B to three Euler angles q = (9, ip, 0) and 
their canonical momenta p = {pe.V^iV^)- In the isotropic case the system reduces further 
to a four- dimensional integrable canonical Hamiltonian system. We follow [ ] but give 
details here to show the system is canonical. As it happens, the transformation is only 
canonical subject to a non-alignment condition. The aligned case is also of interest and is 
treated in subsection 4.2. A reduction similar to the one in the present section was carried 
out in [23] for the nonintegrable anisotropic free rod. 

Let 

' cos 9 cos cos ip — sin sin ip cos 9 cos sin ip + cos ip sin — sin 9 cos ' 
R — — cos 9 sin cos ip — cos sin ip — cos 9 sin <p sin ip + cos (p cos ip sin 9 sin <p 

^ sin 9 cos ip sin 9 sin ip cos 9 j 

be a parametrisation of the rotation matrix R in (2.1) in terms of Euler angles. Here we 
have followed the convention used, e.g., by Love [12], thus 9 is the angle the tangent to 
the rod makes with the magnetic field, ip is the azimuthal angle about the field direction 
and is the twist angle about the centreline of the rod. We may assume, without loss of 
generality, that B is directed along the e% vector of the fixed coordinate system. For the 
triple B we then have 

/ — sin 9 cos d> \ 



B(<Z) = VC 3 R(q)k= ^C 3 



sin 9 sin 
y cos 9 J 



(4.1) 



where k = (0,0, 1) T . On inserting the Euler angles into the strains (2.3) and using the 
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/ 


mi 






I 


Ki{ 


9' sin ( 


f> — ip' sin 9 cos (ft 


)\ 






m = 




m 2 








K 2 ( 


9' cos 


(j) + ip' sin 9 sin 


) 


= Lp, 


(4.2 




V 


m 3 


J 




\ 




K 3 (q 


/ + ^'cos0) 


) 







where 



/ sin 9 sin 6 — cos < 



sin 9 



sin 9 cos ( 




sin ( 




cos 6 1 cos (j) 
— cos 6* sin . 
sin 9 



and p=(p e ,p^, 



J 



are the canonical momenta defined by pg = dW(q,q')/d9', p^ = dW(q,q')/d^', p^ = 
dW(q,q')/d(p', with W(q, q') = W(u(q, q')) in terms of the strain energy function W 
defined in (2.4). Finally, for the force we may write 

n(q,p) = R(q)v(q,p) , 

for some non-constant triple v. Decomposing v as 

V (g, p) = V± (g, p) i ± + W || (g, p) 

where zy and ij_ are unit triples parallel and perpendicular to k, respectively, we obtain 
C 2 = B • n = y^Rk ■ Rv = s/C~ z v ■ (R T R) k = y/C^v\\. (4.3) 

Furthermore, 



RvRv+ s/CzLp .Rk = ^ r + \vl+ ^/C~ 3 p ■ (LR T ) k, 

2C-~ 



which allows us to solve for v±: 
v±(q,p) = 



2C l - C ^-2*jQ i p-{LRT)k 



20,-^-2,/ C 



3 P<p, 



(4.4) 



(4.5) 



where, without loss of generality, we have taken the positive solution. If the vector 
perpendicular to k is taken to be i± = (1, 0, 0) T we obtain 



( 



- sin 9 cos ( 
sin 9 sin 
cos 9 



\ 



C? 



' cos 9 cos (p cos ip — sin <f> sin ip ^ 
— cos 9 sin <p cos ip — cos sin ip 
sin # cos ip 



(4.6) 



c 2 - > 

Note that this transformation is well defined as 2C 1 — — 2vC 3 p^ = > 0. 

Eqs (4.1), (4.2) and (4.6) give the sought transformation. The Jacobian matrix of its 
inverse is 

d(q,p) 



G 



d (m, n, B) 



(4.7) 
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In order for the bracket (3.26) to be transformed to canonical form we need to verify that 
GJG — </, 

where J is the structure matrix defined in (3.21) and J is the standard canonical structure 
matrix in M 6 . We show in the Appendix that this is indeed the case, provided that v± > 0, 

i.e., provided that n and B are not aligned. Without this condition the necessary inverse 

c 2 

transformation does not exist. Note that n and B are aligned if and only if 2C X —-^ = 2m- B. 
Now from (3.26), (3.24) and the conservation of C\ we have 

2^(m-B) = -^(n-n) = 2d 3 -(Bxn), (4.8) 

which vanishes if n and B are aligned. Thus the alignment condition is well defined: if the 
force and the magnetic field are aligned anywhere they are aligned everywhere along the 
rod. 

The Hamiltonian (3.22) is transformed to 

H = ors rl ■ 2 a be sin 2 6 (K 2 - (K 2 - K x ) cos 2 0) + p% {K x sin 2 + K 2 cos 2 0) 
2KiK 2 sin 9 r 

+ p\ (cos 2 9 (Ki sin 2 + K 2 cos 2 0) + (K X K 2 /K 3 ) sin 2 6) 
+ 2{K 2 — KijpQPf/, sin 9 cos 9 sin cos — 2{K 2 — Ki)psp^ sin 9 sin cos(0.9) 
—Ip^p^ cos 9 [K\ sin 2 + K 2 cos 2 



C C 2 
-^cos9 + sinfl cos ip\ 2C 1 - 
'C 3 V C 3 



+ -^cos^ + sin^cos^W2C 1 - -±-2>JC 3 p r (4.10) 



Note that the effect of the magnetic field, encoded by C 3 , has been transferred from the 
Lie-Poisson bracket to the Hamiltonian. The limit C3 — > is singular, as was also observed 
in [14]. 

4-1. The isotropic case 

In the isotropic case (K\ = K 2 = K) the Hamiltonian (4.10) reduces to 

w Pi 1 (^-P^osfl) 2 1 2 C 2 
H= 2K + 2* sin 2 * + 2K^ + 7m C °" 



a 



+ sm9cosipJ2C 1 --±-2 v / C 3 p r (4.11) 



C 3 

Since this Hamiltonian does not depend on the angle the momentum p^ = m 3 is a 
constant. We also have the additional integral (3.30), which in canonical variables reads 

C 

I = ^/C 3 Kcos9 + -j==p^ 




CaP^sin^-cosV ^ Jf^ )) , (4-12) 



rendering the system completely integrable. 
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Hamilton's equations are 



9' = |, (4.13) 

= Qpy, - pt cos 6) _ \fC~z cos ip sin 8 

K sin 2 6 ~ j 2Ci _Z_ 2 ^ 

, (Pi/, cos 9 - p J (p^ - p cos 0) C 2 

Va = ~ 3"^ === Sin 



cos^cos^^/2C 1 -^-2 v ^^, (4.15) 




pi = Bin0sin^W2(7 1 - ^ - 2^/0^. (4.16) 



Helical solutions about B = e 3 have ^ = const., = const. By integrating (4.16) subject 
to these conditions and inserting the result into (4.14) one can show that no helical 
solutions about e% exist in this non-aligned case. 

The phase space of the reduced system (4.13)-(4.16) can be explored by means of 
(planar projections of) two-dimensional Poincare sections of level sets of the integrals 
I and 7i. Fig. 4.1 shows examples in which we have fixed I = 1.00995 and recorded 
intersections of a few orbits with planes of section given by cosi/j = a, for various values 
of a. 



4-2. Alignment of force and field - The superintegrable case 

Here we consider the physically realistic case that the force and the magnetic field are 
aligned. We already showed that if force and field are aligned anywhere then they are 
aligned everywhere. From (3.19) it follows that c£ 3 x B = = d 3 x n, i.e., n is aligned 
with d 3 , while also n and m are constant. This means that solutions are twisted straight 
rods. Hence the aligned case is maximally superintegrable with solutions lying on one-tori. 
Note that this conclusion holds irrespective of whether the rod is isotropic. 



5. Generalised magnetic rods 

5.1. A rod in a nonuniform magnetic field - Four-field model 

By inspection of the structure matrices (3.3), (3.11) and (3.21) we may go to the next 
member of the hierarchy and consider the equations 

m ' + r' xn = 0, n' + r'xB = 0, B' + r' x D = 0, D' = 0. (5.1) 

The equation for B can be integrated to give B x = y, B y = —x, B z = 0, where x, y, z 
and B x ,B y ,B z are components of r and B relative to the fixed frame {ei, 62,63}, and 
we have chosen e 3 in the direction of D. Thus eqs (5.1) can be thought of as describing a 
rod in a linearly-varying magnetic field generated by a uniform 'hypermagnetic' field D. 
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Figure 1. Poincare plots for (4.13)-(4.16) with various sections. In each diagram orbits 
are displayed for energy levels Tt = 1.90, 1.50 and 1.37, while X = 1.00995, A = 0.01, 
Ci = 1.02, C 2 = C 3 = P4, = 1 and K 3 /K = 3/4. 



In the director frame the equations take the Hamiltonian form 
/ m V 



n 
B 

V D / 

and structure matrix 



1 



J (m, n, B, D) VH (m, n) , with H (m, n) = -u ■ m + d 3 ■ n 



J 



Or we can write 

m' 
n' 
B' 
D' 



-J 1 



( m n B D \ 

n B D 

B D 

\ D J 



{ m ^}(m,n,B,D) = m X U + n X d 3 , 
{ n ^}(m,n,B,D) = n X U + B X d 3 , 
{ B ^}( m ,n,B,D) = B X U + D >< d 3, 
{D,^} (m ,n,B,D) = D X U ' 



(5.2) 



(5.3) 
(5.4) 
(5.5) 
(5.6) 
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where the Lie-Poisson bracket is constructed from (3.26) through the addition of another 
semidirect and non-semidirect extension: 

{/> #}(m,n,B,D) = ~ m ' ( V m/ X V m #) - n • (V m / X V n £ + V n / X V m #) 

- B • (V m / x V B g + V B / x V m g) - B • (V n / x V n g) 
-D'(V m /xV D j + V D /xV m5 )- D-(V B /xV nfl ) . 

S V ' V v ' 

evolution of hyperfield effect of hyperfield 

This twelve- dimensional system has four independent Casimirs: 

Ci = m • D + n • B, (5.7) 

C 2 = I B .B + n-D, (5.8) 

C 3 = B ■ D, (5.9) 

C 4 = D ■ D. (5.10) 

In the isotropic case there are now three independent first integrals besides the 
Hamiltonian, 

h = Km-d z (if K 1 = K 2 = K\ (5.11) 

I 2 = n • m + lfB -d 3 (if Ki = K 2 = K) , (5.12) 

/ 3 = in-n + m- B + KD ■ d 3 (if K x = K 2 = K), (5.13) 

making the system completely integrable. If C4 = then D = and the system reduces 
to that of the magnetic rod in the previous section. The system loses rank as the Casimir 
C4 = necessarily implies C3 = and the two Casimirs lose their independent meaning. 
Interestingly, the integral J3 then becomes a Casimir (cf. (3.26)), whose preservation 
does not rely on isotropy anymore. By using the four Casimirs (5.7)-(5.10) the twelve- 
dimensional system can in principle be reduced to an eight-dimensional canonical system. 

Alignment again defines a special case. It can be shown in the same way as in the 
previous section that if D and B are aligned anywhere then they are aligned everywhere. 
From (5.1) it then follows that D and B are aligned with d 3 , which is therefore constant. 
Thus all solutions are twisted straight rods. Again we find that alignment leads to a 
maximally superintegrable case with solutions existing on one-tori. 

5.2. A Lax pair for the isotropic case 

The previous sections show a hierarchy of rod models based on the form of the structure 
matrices. Conditions on the constitutive relations determine if the model is integrable. 
Here we give a compact Lax pair formulation of the integrable hierarchy of isotropic rod 
problems. Consider the parametrised Lax pair 

d r(/i) = [rOi),d 3 /i + ul , (5.14) 



ds 
where 



r fa) = Kd 3 n + r + r^- 1 + . . . + r n/ i" n e so (3) , n e N, 
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with 





/ 





-1 









/ 





-Us 


u 2 


\ 


d 3 = 




i 










and u = 




u 3 





-m 






V 











) 




V 


-u 2 


Ui 





/ 



This Lax pair was proposed in [24] to study monodromy of a generalised family of 
symmetric Lagrange tops. It describes our hierarchy of rod models if we associate the 
terms in the expansion of T by /i with our field variables: r = rh, Ti = n, T 2 = B, 
T 3 = D, etc. The non-canonical equations for the force-free rod (n = 0), Kirchhoff rod 
(n = 1), rod in uniform magnetic field (n = 2) and rod in nonuniform magnetic field 
(n = 3) are obtained by equating like powers of /i in (5.14), while the first integrals are 
generated by 

li = — -residue M=0 (// _1 trace [r (/i) 2 ]) , for i — —1,0, 1, ... ,n — 1, 
the Casimirs by 

Cj = — -residue M= o (// _1 trace [r (/i) 2 ] ) , for % = n, n + 1, n + 2, . . . , 2n, 
and the Hamiltonian is given by 

Zo k-k 3 //_A 2 

K 2KK 3 \K J ' 

6. Conclusion 

We have shown the equations for a conducting rod in a uniform magnetic field to sit 
(as third member) in a hierarchy of rod models, the first two of which are the classical 
force-free rod and the Kirchhoff rod. The hierarchy can be interpreted as describing a 
succession of generalised body forces that are all of the form F x d 3 for some field F, 
i.e., they act normal to the rod. In the case of a transversely isotropic unshearable and 
inextensible rod the equilibrium equations are found to be completely integrable. The 
integrable hierarchy of isotropic rod models is generated by a parametrised Lax pair. The 
integrable models also have maximally superintegrable cases which allow solutions to be 
derived by purely algebraic means, giving straight twisted configurations for the rod. 

The uniform magnetic rod gives a physical realisation to the abstract 'twisted top' 
equations introduced in [14]. The fourth member of the hierarchy also has a physical 
interpretation as a conducting rod in a non-uniform linearly-varying magnetic field. 
The next member is not so easy to interpret as the magnetic field B depends on the 
configuration of the rod; it is no longer an external field. Not all the conserved quantities 
have a physical interpretation, for example (3.26) and (3.30) in the three-field model. This 
is not unusual: the Kovalevskaya integral (3.18) for the heavy top (or Kirchhoff rod), for 
instance, does not have a physical interpretation. 

Unlike the Lagrange case the Kovalevskaya integrable case is destroyed by the 
perturbation due to the magnetic field, although a Kovalevskaya type case may exist 



Integrability of a conducting elastic rod in a magnetic field 



16 



at a different condition. A hierarchy of generalised Kovalevskaya tops is given in [25], 
where also a parametrised Lax pair for this top is derived. Here the perturbation is a 
body moment rather than a body force. The Lax pair was constructed by exploiting the 
Lie-Poisson structure on so(3,n), created by the semi-direct sum of so (3) and n copies 
of R 3 . 

It is interesting to note that Casimirs are 'promoted' to first integrals (and thus become 
conditional on isotropy) as a new field is added in going to the next level of the hierarchy. 
For instance, at the second level of the hierarchy n is added as a uniform field and hence 
|n • n is a Casimir. In the next perturbation, by the field B, the Casimir is perturbed to 
~n • n + m • B. After one more perturbation, by the field D, this Casimir is turned into the 
first integral |n • n + m • B + KD • d 3 . By contrast, the Casimir n • m at the second level 
is perturbed directly into the integral n ■ m + KB ■ d 3 at the next level and remains the 
same one level up. 

We have not discussed boundary conditions of the equations. While rigid-body 
problems are initial-value problems, the rod problems considered here are typically two- 
point boundary- value problems. Buckling results for various types of boundary conditions 
can be found in Wolfe's work [4, 5, 6] and in [11]. It is natural to wonder how other 
perturbations such as extensibility and shearability affect the hierarchy of rod models 
described here. We intend to take up this question in future work. 



Appendix A. The canonical transformation matrix 



Expressing the new canonical variables by inverting the non-canonical variables (4.1), 
(4.2) and (4.6) yields 

,-1 B 3 







tan 



P<t> = m 3 , 



cos , 

m 1 B 2 - m 2 Bi 

B 2 



-1 



B l 



ip = cos 



-1 



n 3 - C 2 B 3 /C 3 



y/(l - B'i/C 3 ) (2Ci - Cf/C 3 - 2 (m 1 B 1 + m 2 B 2 + m 3 B 3 )) 
Pip = {miBi + m 2 B 2 + m 3 B 3 ) / y/C^. 
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The transformation matrix G in (4.7) is given by 



G 













































B'i+Bl 


Bi+Bi 


#31 


#32 


#33 


934 


935 


936 


937 


938 


B 2 
















rn 2 


ra\ 






^jGz-Bl 










1 


















c. 



1 — Co Bo 



\ Bjy/Vz B 2 /VC~3 B 3 /VC~ 3 Tru/VUs" m 2 /^C~ 3 





#39 
{Cs-Blf 2 



m 3 /^/Cl 



where 



#31 
#32 

#33 
#34 
#35 

#36 

#37 

#38 
#39 



Bi {n 3 


— C 2 B 3 /C 3 ) 




A 


B 2 (n 3 


— C 2 B 3 /C 3 ) 




A 


S 3 (n 3 


— C 2 B 3 /C 3 ) 




A 


0. 




0, 




1 








mi (n 3 


— C 2 B 3 / C 3 ) 




A 


m 2 (n 3 


— C 2 B 3 / C 3 ) 



A 

C 2 (2^ - C 2 2 /C 3 - 2m ■ B) /C 3 
A 

(n 3 - C 2 S 3 /C 3 ) (S 3 (2C*! - C 2 2 /C 3 — 2m • B) — m 3 (1 - B 3 /C 2 )) 



A 



with the denominator given by 



A = (2C 1 - Ci/C 3 - 2m • B) y (1 - B 2 3 /C 3 ) {2C 1 - Cl/C 3 - 2m • B) 
The non-canonical structure matrix is given by 
/ 



C 2 B 3 /C 3 ) 2 . 



J 






-m 3 


m 2 





-n 3 


n 2 





-s 3 


s 2 


\ 


m 3 





—mi 


n 3 





-m 


s 3 





-Bi 




-m 2 


nil 





-n 2 







-s 2 


Sx 










-n 3 


7? 2 





-s 3 


s 2 













n 3 





-m 


s 3 





-Sx 













-n 2 


ni 





-s 2 


S x 



















-S3 


s 2 






















B 3 





-S x 






















-B 2 


Si 























/ 
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while the canonical structure matrix is given by 



J = 



( 











1 





\ 
















1 






















1 




-1 






















-1 














V 








-1 
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